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ABSTRACT 

We study the possibility of quasar outflows in clusters and groups of galaxies heating 
the intracluster gas in order to explain the recent observation of excess entropy in this 
gas. We use the extended Press-Schechter formalism to estimate the number of quasars 
that become members of a group of cluster of a given mass and formation epoch. We 
also estimate the fraction of mechanical energy in the outflows that is imparted to 
the surrounding medium as a function of the density and temperature of this gas. We 
finally calculate the total amount of non-gravitational heating from such outflows as 
a function of the cluster potential and formation epoch. We show that outflows from 
broad absorption line (BAL) and radio loud quasars can provide the required amount 
of heating of the intracluster gas. We find that in this scenario most of the heating 
takes place at z ~ 1-4, and that this "preheating" epoch is at lower redshift for lower 
mass clusters. 

Key words: Cosmology: Theory — Galaxies: Intergalactic Medium — X-rays: Galax- 
ies: Clusters 



1 INTRODUCTION 



Clusters and groups of galaxies contain a large amount of 
hot gas, besides galaxies and the gravitationally dominant 
dark matter. This hot X-ray-emitting gas known as the in- 
tracluster medium (ICM) represents a part of the baryonic 
matter of the universe that is not associated with individual 
galaxies but remains trapped in the deeper gravitational po- 
tential of galaxy clusters. Hierarchical models of structure 
formation have been very successful in explaining many ob- 
served properties of galaxies and galaxy clusters. Neverthe- 
less, some puzzling problems remain open and unexplained. 
Models of cluster formation in which the intergalactic gas 
simply falls into the dark matter dominated gravitational 
potential well (so-called infall models) fail to reproduce all 
the structural properties of the local cluster population (e.g., 
Evrard & Henry 1991; Navarro, Frenk & White 1995; Mohr 
& Evrard 1997; Bryan & Norman 1998). There is certainly 
some additional physics driving ICM evolution. 

Recent X-ray observations have provided evidences for 
some non-gravitational heating of the diffuse, high density 
baryons in the potential wells of groups and clusters of galax- 
ies, in addition to the heating during the gravitational col- 
lapse. One of the first evidences was in the shape of the 
L x — T relation, which is steeper than the self-similar be- 
haviour L x oc T 2 predicted in the case of gravitational 
processes only. As early as the emergence of ROSAT and 
Einstein data, several authors proposed that the missing 
element is the existence of a "preheated high-entropy" in- 



tergalactic gas prior to a cluster's collapse (David et al. 1991; 
Evrard & Henry 1991; Kaiser 1991; White 1991). Later, Pon- 
man et al. (1999), and Lloyd-Davies et al. (2000) found direct 
evidence of an entropy excess with respect to the level ex- 
pected from gravitational heating in the centre of groups. 
Ignoring the constant and logarithms, one can define the 
"entropy" as S = T/n 2 ^ 3 . The excess entropy, or equiva- 
lently, the excess specific energy, flattens the density pro- 
file decreasing the X-ray luminosity which is proportional to 
the square of the density. The effect is stronger in poorer 
clusters, where the excess energy associated with the excess 
entropy is comparable to the gravitational binding energy, 
while rich clusters, where gravity is dominant, are mostly 
unaffected. This produces a steepening of the L x — T rela- 
tion. 

The most popular scenario to successfully explain these 
thermal properties of the ICM has been the "preheating" 
scenario. For this scenario, the candidate processes which 
have been looked into are strong galactic winds driven by 
supernovae. However Valgeas & Silk (1999) showed that the 
energy provided by supernovae cannot raise the entropy of 
intergalactic medium (ICM) up to the level required by cur- 
rent observations. The observed amount of required energy 
injection depends on the epoch, and have been in the range 
of 0.4 - 3 KeV per gas particle (Navarro et al. 1995; Cava- 
liere et al. 1997; Balogh et al.1999; Wu et al. 2000; Lloyd- 
Davies et al. 2000, Borgani et al. 2001). For example, Wu et 
al. (2000) showed that galactic winds can impart only < 0.1 
keV per particle. Moreover, Kravtsov & Yepes (2000) es- 
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timated the energy provided by supernovae from the ob- 
served metal abundance of ICM and found that the heating 
only by supernovae driven outflows requires unrealistically 
high efficiency. On the other hand, quasar outflows may be 
much more powerful and plausible candidates of the heating 
(Valageas & Silk 1999). In this paper, we focus on the role 
of quasar outflows in this regard. 

The epoch of the energy input also remains uncertain. 
Lloyd-Davies et al. (2000) put an upper limit of z max ~ 7-10 
on the preheating epoch, from their estimate of excess en- 
tropy in groups. For AGNs, there have been no additional 
constraints like the metal abundance in the case of super- 
nova heating (Kaiser & Alexander 1999). Recently, Yamada 
& Fujita (2001) have looked into the deformation of the cos- 
mic microwave background (the Sunyaev-Zel'dovich effect) 
by hot electrons produced at the shocks produced by jets 
from AGNs. They showed that the observed excess entropy 
of ICM and COBE/FIRAS upper limit for Compton y- 
parameter are compatible with each other only when the 
heating by jets occurred at relatively small redshift (z < 3). 
Thus they questioned the "preheating" scenario as their re- 
sult suggests that the heating occurred simultaneously with 
or after cluster formation. 

In this paper, we calculate the heat input from quasar 
outflows inside clusters. We calculate the mechanical work 
done by various kinds of quasar outflows and the excess en- 
ergy imparted by them onto the intracluster medium via 
pdV work. For the statistics of quasars inside clusters, we 
use the extended Press-Schechter formalism. Finally we cal- 
culate the excess energy per particle and tally them with 
available observations. 

In the next section, we discuss the abundance of quasars 
inside clusters of a given mass. We then discuss the evolution 
of quasar outflows and the mechanical work done by them 
in §3. We use these concepts to calculate the heating of the 
ICM in §4. We then discuss the implications of our results 
in §5. Throughout the paper we assume a flat universe with 
a cosmological constant, with Qo = 0.35, Oao = 0.65, and 
h = 0.65. 



2 QUASARS INSIDE CLUSTERS 

For a proper evaluation of the heat input from quasars inside 
clusters, one first needs to calculate their abundance and its 
dependence on the quasar mass, cluster mass, and the clus- 
ter formation redshift. Observationally, it is still difficult to 
obtain good statistics of quasars inside clusters. Estimation 
of the galaxy-QSO correlation function have shown that at 
low redshifts (z Ss 0.4), quasars typically reside in small to 
moderate groups of galaxies and not in rich galaxies (e.g., 
Bahcall & Chokshi 1991; Fisher et al. 1996). It is still un- 
certain whether or not their optical or radio luminosities 
depend on their environments. Some studies (e.g., Elling- 
son, Yee & Green 1991) have shown that radio-loud quasars 
preferentially reside in richer clusters at higher redshift, al- 
though more recent studies (e.g., Wold et al. 2001) do not 
find any such correlation. 

There has been, however, considerable work in relating 
the observed quasar luminosity function or the radio lumi- 
nosity function (for radio-loud quasars) with the mass func- 
tion of galaxies as prescribed by the Press-Schechter (PS) 



formalism (e.g., Haehnelt & Rees 1993; Haiman & Loeb 
1998; Yamada, Sugiyama & Silk 1999). 

At z = 0, Yamada, Sugiyama & Silk found that one can 
reproduce the abundance of the radio sources powered by 
AGNs (i.e., leaving aside the radio sources powered by star- 
bursts) by assuming that a fraction f r of the halos from PS 
formalism become radio-loud quasars, where, f r ~ 0.01 for 
M h ;> 10 12 M Q , and f r = for M < 10 12 M Q . They assume 
an upper limit on Mh of 10 14 Mq. Since it is known that 
radio loud quasars constitute a fraction 0.1 of the quasar 
population (Stern et al. 2000) , this means that a fraction 
f q = 10/r of the halos from PS formalism become quasars. 
In other words, /, ~ 0.1 for 10 12 < M h < 10 14 M , and 
f q = 0, otherwise. Yamada, Sugiyama & Silk (1999) assumed 
this fraction to be a constant for all redshifts. In this model, 
the rate of formation of quasars is given by the derivative of 
the PS mass function at the relevant mass scale. 

This is similar to the model adopted by Haiman & Loeb 
(1998) and Furlanetto & Loeb (2001; FL01). FL01 showed 
that at high redshift (z <; 4) the rate of formation of quasars 
is a fraction f q ~ 0.1 of the rate of formation of halos from 
PS formalism, if a life time of order 10 7 yr is assumed for 
the quasar. There is, however, a difference, in that they had 
only a lower limit to Mh and no upper limits (see eqns (2) 
& (3) of Haiman & Loeb 1998). They mention that at low 
redshifts their formalism does not predict any decline as ob- 
served in reality, and for which they consider their model 
only at high redshift. It is possible that this mismatch is 
due to the lack of upper limits, since the the differentia- 
tion of PS mass function for objects with an upper limit 
in mass decreases at low redshift, (Haiman, Z. 2001, private 
communications). In any case, Haiman & Loeb (1998) found 
that this prescription yields a matching quasar luminosity 
function that is observed, at redshift z > 2.5. We find later 
that most of the heating of the ICM gas (even for our least 
massive cluster) occurs at z > 2 (Figure 7). We will there- 
fore assume for simplicity that this fraction f q ~ 0.1 at all 
redshifts (as in Yamada et al. 1999). 

In this paper, we would like to have a conservative esti- 
mate of the quasar abundance. Also, we would like to calcu- 
late the abundance of quasars in clusters including low mass 
groups of galaxies. For this reason, we assume the value of 
f q as above, but use an upper limit of 10 13 Mq. Since the 
mass function decreases steeply at the higher mass end, this 
should not change the value of f q substantially. In brief, we 
assume that, 

/ 0.1 if 10 13 k, M h Z 10 12 M Q 
h ~ \ if M h > 10 13 Mq , Mh < 10 12 M Q , 1 ' 

motivated by the model of Yamada, Sugiyama & Silk (1999), 
and by the fact that a similar prescription by Haiman & 
Loeb (1998) recovers the quasar population at high redshift, 
and relate the rate of formation of quasars with that of halos 
in the PS formalism. 

We are, however, concerned with the statistics of 
quasars inside clusters. For this one needs to have an ex- 
tension of the PS mass function which can predict the prob- 
ability of a given halo becoming a part of bigger object later, 
or the probability of an object having had a progenitor of 
a given mass at an earlier epoch. Such extensions of the 
PS theory have been studied in detail by Bower (1991) and 
Lacey & Cole (1993), for example. 
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In the standard PS theory, the mass function, i.e., the 
fraction of regions with mass in the range M, M + dM which 
have overdensity <5 in excess of S c (which is the threshold for 
perturbations becoming non-linear), is given by, 



f PS =f(M,S c )dM = 



2^ «) 3/2 



exp 



2< 



dM 



dM.(2) 



Here, ctm is the mass variance of the perturbation at the 
mass scale M. The relation between the number density of 
objects in the mass range Al, M + dM with f(M)dM is, 

n(M)dM = ^f(M)dM , (3) 

where p denotes the background mass density. 

In the extended PS theory, the fraction of regions of 
mass M, contained within a larger scale region of mass M' 
and overdensity 5' , which are more overdense than S c , is 
given by, 



f(M, 5 c \M',6')dM 



(8 C - 5') 



2tt (er 2 



x exp 



M 

(S, 



- a 2 ,W 2 



8'f 



2(<r 



dM 



dM .(4) 



This expression recovers the simple PS mass function in the 
limit M' — ► oo and 8' — > 0, relevant for the whole universe. 

If we then identify M' with M c i , the mass of a cluster, 
and 8' = 8 c (zf), the threshold overdensity of the cluster at 
its formation epoch z/, we can then obtain the mass fraction 
of M c i which have been parts of progenitors of a given mass 
range (M to M + dM) at a given (earlier) redshift. If we also 
identify this mass range with that of the quasars as in the 
standard PS theory (10 12 -10 13 M Q ), and use the fraction of 
these halos that become quasars (f q ; eqn ^), we will obtain 
the mass fraction of the final cluster which have been quasars 
at some given earlier epoch. 

To obtain the rate of formation of these quasars (inside 
a future cluster of mass M c i), we should differentiate the 
above expression. A simple differentiation will, however, not 
give the correct result, since there will be a negative con- 
tribution from the merging of halos out of this mass range. 
One would get a negative rate of formation of such quasars 
at some point if the rate at which they disappear beyond 
this mass limit is not taken into account in a proper man- 
ner. This 

Consider the abundance of objects in a given mass range 
M,M + dM at two successive epochs zi and 22 {zi < 
The abundance f(M)dM at Z2 will be given by /2 = /1 + 
F — D, where /2 and /1 are the abundances at epochs Z2 
and zi, the term F denotes the abundance of newly formed 
objects in this mass range during the epoch Zi and Z2 (from 
merger of smaller objects), and D signifies the abundance 
of objects that moved out of this mass range as a result of 
merger (into bigger objects). A simple differentiation of the 
PS function, involving the difference (/2 — fi) will therefore 
depend on both F and D. For a given range of mass, D 
is very small at a very early epoch, but it increases with 
time (see, e.g., Fig. 5 of Haiman & Menou 2000), and at a 
later epoch can become larger than F. Therefore, at lower 
redshifts, a simple differentiation can imply a negative rate 
of change of abundance. If the contribution of D is neglected, 
one would then incorrectly get a negative value of F. 

This problem has been encountered in the case of or- 



dinary PS function by many authors. While studying the 
rate of mergers in the context of background radiation from 
starbursts and AGNs, Blain & Longair (1993) noted that a 
simple differentiation of the PS function leads to a negative 
rate of formation of objects in a given mass range. In other 
words, the actual rate of formation of objects is given by, 



fform — /PS + frr 



(5) 



They performed a simulation assuming a simple power spec- 
trum and obtained a fit for the rate of these objects merging 
to form bigger objects. They found that fmerger can be ap- 
proximated well (in the Einstein de-Sitter universe) by, 



fmerger — , exp 



(6) 



'2<7 2 _ 

where the value of a ~ 1.35 and <f) ~ 1.3—1.7. This problem 
has also been investigated by Sasaki (1994) and Percival 
& Miller (1999). With the extension of the PS formalism, 
one can now calculate this merging rate (see also, Chiu & 
Ostriker 2000). 

We note here that in the case of the merger of a lower 
mass quasar into a more massive quasar, our implicit as- 
sumption is that the central black holes also merge and form 
a bigger black hole appropriate for the bigger quasar (see 
below; eqn|l^). In other words, we assume that the central 
black hole mass always traces the halo mass. 

The rate at which an object of mass M at epoch z 
merges to form a bigger object of mass M' is given by (Bower 
1991; Lacey & Cole 1993), 



d 2 



dM'dz 



(M -> M'\z)dM' = 



1 



2tt o 2 M ,{a 2 u 



\3/2 



x exp 



S c (z) 2 (al 



dM 



I d5c(z) 
dz 



\dM . 



Here, p(M — » M'\z)dM' is the probability of an object of 
mass M merging to become an object of mass within the 
range M' , M' + dM' at redshift z. The rate of disappearance 
of objects of a given mass, hmerger, should be essentially, 
fif r°° j2„ 

Ll J merger 



dz 



-(M, 



fps(M,z) 



ri 2 . 
dM'dz 



(M-> M'\z)dM' ,{8) 



where the lower limit of the integration is chosen to be such 
that the merged object is at least twice as massive as the 
merging object. We show the result of this integration for a 
sCDM universe (with a COBE normalized power spectrum) 
in Figure 1, and show the fit of Blain & Longair with <f> — 0.9 
and a — 1.35. We find that the merging rate is fit by a lower 
value of <f) than they assumed, although the difference is a 
factor of order unity. It is possible that this difference is 
due to the specific assumption in the simulation done by 
Blain & Longair (1993), e.g., in the power spectrum being a 
simple power law (Blain, A. 2001, private communication), 
or it can be a result of the lower limit (2M) chosen by us. 
At any rate, if we chose the above integral to represent the 
merger rate then it would be a conservative estimate, since 
decreasing the lower limit would simple increase the value of 
the merger rate, and in turn, the formation rate of objects. 

We have found that the result in the case of the ACDM 
universe can be fit by a similar function, with <f> being re- 
placed by 0.9 —j^-, for a. = 1.35 with an accuracy of order 
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Figure 1. The rate of disappearance of objects is compared 
in different cosmologies. The thin solid line shows the result of 
integration in eqn (H) for a SCDM universe for M = 10 14 M© 
and the dotted line shows the Blain &: Longair (1993) fit with 
<j> = 0.9 and a = 1.35. The thick solid line shows the result of the 
integration in the ACDM universe and the dashed line refers to a 
fit described in the text. Both curves use 4-year COBE normalized 
spectra. 



Sa 5%. We, however, do not use these fits in our calculation, 
and evaluate the integral numerically for our purpose. 

To be precise, this rate of disappearance is valid for 
the objects following the PS mass function, i.e., for objects 
which are not already parts of bigger objects. Motivated by 
the extension of the PS formalism, we here posit that the 
rate of disappearance of objects inside a bigger object also 
has the same form, with fps(M, S c ) in eqn (H) being replaced 
by f(M,8 c \M',S'). There is admittedly no way of verifying 
the truth of this ansatz at present, since this would involve 
more extensions of the PS theory than that is available now. 
It will also involve comparing the merger rates inside and 
outside of clusters. It, however, leads to a conservative es- 
timate for the formation rate of quasar in a cluster. As the 
work of Bower (1991) has shown, growth of perturbations 
inside a cluster is enhanced compared to in the field. This 
means that the merging rate of objects of a given mass in- 
side a cluster should be larger than that in the field. Here, by 
assuming a comparable merging probability (the factor that 
multiplies the abundance of objects, given by fps), we are in 
a way underestimating the rate of disappearance (f merger), 
and in turn, the rate of formation (// orm ) of quasars in a 
cluster. The final result of total heat input from our formal- 
ism should, therefore, be a conservative estimate. 

We show the results of adding the rate of disappear- 



ance in Figure 2. The dashed lines, which show the term 
— {dfps/dz) (or equivalenlty, d ^ s ), become negative at 
lower redshift, suggesting the need for the addition of the 
rate of disappearance of objects. The dotted lines show the 
result of adding this rate, using the Blain & Longair (1993) 
fit (with 4> = 0.9), and the solid lines show the result using 
the integral in equation |j| The upper panels of the figure 
show the case of objects (of mass 10 13 M©) in the field, and 
the bottom panels show the case for these objects inside a 
cluster of Mel = 10 Mq (z/ = 0). For these curves, we have 
used the integration in equation g, and the extension of the 
PS mass function for the abundance of objects inside the 
cluster, as explained above. The left and right panels show 
the cases for sCDM and A-CDM universe. As the bottom 
panels show, the addition of the disappearance rate does 
not suffice to make [— (dfps/dz)] (inside clusters) positive 
at low redshift. This is suggestive of the enhanced growth of 
perturbation, and the need for a larger rate of merger inside 
clusters, as previous authors have noted. 

We also show the case for M c i — 10 14 M© by the long 
dashed line in the lower right panel of Figure 2. Comparison 
with the solid line (for M c i = 10 15 M©) shows that the for- 
mation rate of these galaxies inside a lower massive cluster 
is larger. This follows simply from the extension of the PS 
formalism (from the dependence on the term (cr^ — It 
is interesting to note that this is consistent with the sugges- 
tion from observation that quasars are preferentially located 
in groups of galaxies instead of rich clusters (e.g., Bahcall & 
Chokshi 1991; Fisher et al. 1996). 

If our formalism is used without any correction, this 
will lead to subtraction of energy input in the final result. 
We circumvent this problem by putting [—(dfps/dz)] = 
when this term turns negative. This will, therefore, provide 
a lower limit to the total energy input from quasar outflows 
in a cluster. 

We can finally write down the rate of formation of 
quasars in a given mass range M, M + dM inside a clus- 
ter of a given mass, M c i, in the form of the rate of increase 
of the fraction of mass of M c i which is in quasars at an epoch 
z, as (remembering the equation^, 



df q , c i(M,z\M c i,z f ) 
dz 



dM = /< 



df(M,S c (z)\M cl ,S c (zf)) 



+f v f(M,5 c (z)\M cl 



dz 

5 c (zf))dM 



dM 



x 



dM'dz 



(M — > M'\z)dM' , (9) 



with the condition that 

dn n 



dn,j 
dz 



where 



= 0, for z < 

-j7r\z<z rl > 0. The integral on the right hand side is eval- 
uated using eqn(Q). Here we have also changed the upper 
limit of the integration to M c ;. The integrand is a rapidly 
decreasing function of M and the value of the integral de- 
pends mostly on the lower limit. 

Here, the threshold density contrast in a cosmological 
constant dominated universe is given by a fit given by Ki- 
tayama & Suto (1996), 



8 c (z) = 1.68[g(z = 0)/g(z)][l + 0.0123 log fi m (z)] . 



(10) 



In our calculations, we have used a fit for g(z) from Carroll, 
Press & Turner (1992), 



,fi A (z)) 



© 0000 RAS, MNRAS 000, 000-000 



Heating of the intracluster medium by quasar outflows 5 



0.1 

0.08 
^ 0.06 

^ 0.04 

T 

0.02 


0.06 

7 0.04 

n 

\ 

3 

o 

T3 

T 0.02 







- 

- 


sCDM 1 

: 


- 


/,!,,,!,,, 


2 


4 6 f 




(1+z) 


I 

- 


I i i i I i i 

sCDM, M cl - 10 15 Mg 




.' / / ^ n\ 
// / V 

// / 
•/ 1 




/ / 


i ;'/ 


' / 
/ 

,1/11,1,,, 



4 

(1+z) 




Figure 2. The net formation rate of objects ( ^£2ng ) is plotted for objects of mass 10 13 Mq (for dM = M), 
in clusters (bottom panels, with M ci = 10 Mq) and in general (top panels), for sCDM (left panels) and 
ACDM (right panels) models. Dashed lines show the term d ^f s , dotted lines show — Lsnn using the Blain 



& Longair fit (with rj> = 1.0), and solid lines show 



■Vfo 



line in lower right panel shows the case for M c i = 10 Mq 



using the integral in equation H. The long-dashed 



" 2[n m ( z )^ - n A (*) + (i + ^i)(i + 2aM)] ' 

where (Lahav et al. 1991), 

Q m (z) = Q m0 (l + zf /[Q mQ {l + zf + Q A0 ] , 

q a (z) = n A0 /[n m0 {i + z) 3 + n A0 ]. (12) 

We emphasize here that the above formalism leads to 
a conservative estimate of abundance of quasars in a clus- 
ter (and, therefore, the final heat input), because (a) we 
ignore the increased pace of growth of perturbation and the 
merging rate inside a cluster, and (b) the lower limit of the 
integration could in reality be smaller than 2M, which is 
probably the reason the rate of formation still turns nega- 
tive at low redshifts even after the addition of merger term. 

To summarise the work in this section, we have used 
the existing ideas for relating the quasar formation rate to 
the Press-Schechter mass function, to estimate the rate of 
formation of quasars inside clusters (as a function of cluster 
mass and formation redshift), utilizing the extensions of PS 
formalism. First, we use eqn ([ij) to relate the PS mass func- 
tion to quasar abundance. The standard PS mass function 
(eqn |^) is then replaced by its extension (eqn Q). Further- 
more, we add the contribution due to merger (into larger 
objects) (eqns ^ & |^) , again using the extensions of PS mass 



function. We finally have the rate of formation of quasars 
inside clusters as given by eqn (^). 

We use this formalism to calculate the total energy in- 
put from outflows from quasars in a cluster, or, a lower limit 
to it. We next discuss the energy input from individual out- 
flows, which we will combine with our calculation of forma- 
tion rate of quasars for the final result. 



3 WORK DONE BY QUASAR OUTFLOWS 

In this section we calculate the energy input from quasar 
outflows into the ambient medium. We model the outflows 
as they evolve in the ambient medium and calculate the pdV 
work done by the outflows. To begin with, we discuss differ- 
ent kinds of outflows that we consider and the characteristics 
of the hosts of quasars. 

3.1 Quasar outflows 

We consider two major types of quasar outflows. For radio- 
loud quasars (RLQ), the outflow is in the form of a tightly 
collimated jet, which deposits energetic particles into a co- 
coon which expands against the surrounding medium. These 
outflows are characterised by the kinetic luminosity of the 
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jet, Lfe. According to Willot et al. (1999), this is correlated 
with the bolometric luminosity Lboi of the quasar, and that 
0.05 S5> Lk/Lboi 1.0. We follow FL01 in arguing that since 
Lboi ~ IOLb (Elvis et al. (1994)), the rest-frame B-band 
luminosity, Lk ~ Lb- 

Radio-loud quasars, however, constitute only about 
10% of the total population of quasars (Stern et al. 2000) . 
We therefore define a factor f for the fraction of quasars 
with outflows, and define f ~ 0.1 for our RLQ model. The 
fraction f a here denotes the number of radio relics/lobes per 
halo, since a radio loud quasar may have several outbursts of 
radio activity. This fraction is therefore a very conservative 
estimate since it is obtained from observed radio luminos- 
ity function and does not take into account the existence of 
radio relics in clusters. 

Another important kind of outflows are encountered 
in broad absorption line (BAL) quasars. The absorption 
troughs are thought to be due to absorbing clouds flowing 
out of the quasars with velocities up to 0.1c. Although they 
are encountered in about 10% of the quasars, it is believed 
that all quasars have such outflows (all the time) and the 
covering fraction of the BAL outflows is about 10% (Wey- 
mann et al. 1991; Weymann 1997). Some authors also believe 
that BAL outflows have a limited lifetime (especially the low 
ionization BALs) and that they have a large covering frac- 
tion in the early phase of a quasar (Voit et al. 1993). For our 
calculation, we use a fraction f a ~ 1 for the BAL outflows. 
We discuss the effect of the uncertainty in these factors on 
the final result in §5. 

We model the BAL outflows as having a kinetic lumi- 
nosity Lt- Following FL01, if Nh is the column density 
of the absorbing gas, f c the covering fraction, and Rbal 
is the size of the absorption system, then Lk is related to 
the outflow velocity vbal as Lk ~ 2Tvf c NHm p RBALV% AL . 
The observed range of these parameters are as follows: 
vbal ^ 0.1c, f c ~ 0.1 (Weymann 1997; but see above), 
Rbal ~ l-500pc and N H ~ 10 22 -10 23 cm -2 (Krolik 1999; 
Gallagher et al. 1999). For these values, the magnitude of 
Lk is close to that of Lb- FL01 also argued that for BAL 
winds Lk ~ 0.1-100Lb, and finally assumed Lk ~ Lb - Since 
this estimate depends crucially on a number of uncertain 
parameters (for example, the fact that the absorption col- 
umn density in optical measurement is much smaller than 
the above mentioned X-ray column density), it may not re- 
ally be a conservative estimate, but it does provide a simple 
scaling which we hope is not too unreasonable. The esti- 
mate Lk ~ Lb ~ O.lZ/Edd is probably not a conservative 
estimate, but an upper limit, in that for a covering fraction 
of 10%, the mechanical luminosity of the wind could not 
be larger than 0.1 of the Eddington rate. Keeping all these 
uncertainties in mind, we assume that Lk ~ Lb for BAL 
outflows. 

We then need to connect Lb of a given quasar with the 
properties of its halo. Firstly, as Haiman & Loeb (1998) have 
shown, the mass of the black hole at the centre is related to 

Lb, as, 

Mb " = 0^93 ( 1.4 x 1033 ergs-0 MQ (13) 

assuming that the quasar radiates at the Eddington lumi- 
nosity. The factor of 0.093 reflects the fraction of the Ed- 
dington luminosity radiated in the B-band, taken from the 



median quasar spectrum of Elvis et al. (1993). Statistically 
speaking, we therefore assume that a fraction f q (eqn(^)of 
all black holes radiate at the Eddington rate for a life time 
of ~ 10 7 yr, while the rest does not radiate at all. 

Secondly, the correlation between the central black hole 
mass and the total baryonic mass of the galaxy (Magorrian 
et al. 1998; Gebhardt et al. 2000) gives M B h ~4x 10" 4 M h , 
where Mh is the total mass of the galaxy, using a value 
of Mbh /Mbaryonic ~ 2-3 x 10~ 3 , and M baryonic /M h ~ 

n b /n ~ 0.2. 

As far as the collimation is concerned, the geometry 
of BAL outflows is still uncertain, whereas the radio jets 
are well collimated. Since some models do suggest a modest 
collimation even in BAL outflows, with a covering fraction 
of f c ~ 0.1 (Weymann 1997), we use the idea of collimated 
outflows for outflows from both radio-loud and BAL quasars. 

We next discuss the evolution of the outflows from 
radio-loud quasars and calculate the fraction of its total 
kinetic luminosity that it deposits into the surrounding 
medium in the form of pdV work. For concreteness, we will 
assume that the corresponding fraction for BAL outflows 
has similar values, and keep in mind the uncertainty in the 
geometry and energetics of BAL outflows. 

3.2 Evolution of outflows 

The standard scenario for outflows from radio loud quasars 
involves a 'cocoon' surrounding the core and the jet, and 
consisting of a shocked ambient medium and shocked jet ma- 
terial (Scheuer 1974; Blandford and Rees 1974). Begelman 
& Ciofli (1989) constructed a simple model of the evolution 
of a cocoon in which the cocoon is overpressured against the 
ICM. In their model, the expansion along the jet axis is de- 
termined by the balance of the thrust of the jet and the ram 
pressure, whereas the thermal pressure of the cocoon drives 
along the direction perpendicular to the jet axis. Results of 
numerical simulations agree with this scenario (Loken et al. 
1992; Ciofli & Blondin 1992). 

Here we adopt the model of the evolution of cocoons fol- 
lowing the approach of Bicknell et al. (1997), which is based 
on the Begalman & Ciofli (1989) model but includes the 
pdV work done by the cocoon, in order to find the frac- 
tion of total energy lost by the quasar to the ICM through 
mechanical work (pdV work). Bicknell et al. (1997) derived 
this fraction (f p dv) to be f p dv = 3/8, for a homogeneous 
ambient medium (their equation (2.13)). 

We will, however, calculate this fraction from numerical 
solution of the equations governing the evolution of the co- 
coon, for the following reasons. Firstly, the derivation men- 
tioned above implicitly assumes that the mean pressure av- 
eraged over the hotspot region is equal to the mean lobe 
pressure (in the language of Bicknell et al. (1997), this means 
£ ~ 1). In fact, their equation (2.13) shows that for constant 
C, one has in general, Pc^f = (l + 2£)Lj/8C, which recovers 
the fraction 3/8 for £ = 1. In reality, however, this ratio does 
not remain a constant in time. Secondly, this derivation is 
valid only during the period when the jet is active. Even after 
the jet switches off, the cocoon, however, continues to evolve 
as a result of its overpressure until it reaches an equilibrium 
pressure with the ambient medium (see also Nath 1995). The 
cocoon, therefore, continues to do pdV work even after the 
jet switches off, and the inclusion of this process will lead to 
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an upward revision of the fraction of total energy that is lost 
in pdV work. Besides, there seems to be some confusion in 
the literature regarding the fraction. For example, Inoue & 
Sasaki (2001) have recently adopted a fraction f pdv = 1/4 
in their calculation of energy input into the surrounding gas. 

We therefore calculate this fraction by numerically solv- 
ing the equations of cocoon evolution. 

We consider two collimated steady jets advancing into 
the ambient ICM. The thermalized jet matter and the shock- 
compressed ICM matter form a cocoon around the jets and 
the cocoon expands with shocks advancing in directions both 
parallel and perpendicular to the jet axis. After this stage 
of evolution, when the jet turns off after a lifetime of tu / e ~ 
3 x 10 7 years (Kaiser 2000) , cocoons still retain high pressure. 
They cool radiatively and expand due to its overpressure till 
it reaches a pressure equilibrium with the ambient medium. 
Thus the relevant equations are: 



dr h 
dt 



dr c 
lit 
dE c 
dt 



\A h p a pcJ 

(PcY 12 

I I , t> ttife 

V Pa / 

/ p \ 1/2 

© ■ 



= U 



dV c 
1 dt 



(14) 
(15) 
(16) 



where Lj is the jet luminosity, p a is the density of the am- 
bient medium and (3c is the velocity of the jet material. As 
the jet is highly relativistic, ft ~ 1. The averaged hotspot 
area Ah ~ 30 kpc 2 (Bicknell et al. 1997) is assumed to be 
larger than the radius of the jet, according to the 'dentist's 
drill' model of the jet (Scheuer 1982). Here ru is the length 
of the jet or the distance of the hotspot from the centre of 
the galaxy, r c is the half-width of the cocoon at the centre 
and V c is the volume of the cocoon given by V c = e v (2-Kr\)rh 
where e v is the geometrical factor depending on the shape 
of the cocoon. For our calculations we have taken the shape 
of the cocoon as biconical and so e v ~ 1/3. Finally P c is 
the pressure inside the cocoon given by P c = (E c /V c )(-y — 1) 
where 7 = 4/3 and E c is the total energy inside the cocoon 
given by E c = Ljt till t < t U f e and E c = Ljtuf e afterwards, 
where tuf e is the lifetime of the jet. 

This is admittedly a simplified model of the evolution of 
the cocoon. In reality, after the jet switches off, one expects 
Rayleigh- Taylor and Kelvin-Helmholtz instabilities to dis- 
tort the cocoon, giving rise to 'buoyant' plumes. This phase 
of the evolution of the cocoon, and its effect on the ambi- 
ent medium, has been recently addressed by various authors 
(Gull & Northover 1973; Churazov et al. 2000; Briiggen & 
Kaiser 2000), mainly with the help of numerical simulations. 
It is possible that this phase adds substantial heating to the 
intracluster medium, but short of doing a numerical simu- 
lation it is difficult to assess its importance. We have also 
neglected the loss of energy through radiation, as modeling 
such losses would involve detail knowledge of various param- 
eters (e.g., electron energy spectrum and the possibilities of 
re-acceleration of electrons) . With the uncertainties involved 
in modeling these process, it seems reasonable to adopt the 
above simplified picture as a pointer and keep the uncer- 
tainties in mind while discussing the final result. In light of 




-5 

Log 10 n (particles/cc) 



Figure 3. The dependence of f p dv with the ambient density n is 
shown for various ambient temperatures T for various (erg/s) 
with and without radiation loss. Solid and dashed lines are for 
Lk = 10 46 erg/s and 10 47 erg/s respectively without radiation 
loss. Dotted lines are for L k = 10 46 with radiation loss. 



this discussion, we will also calculate the final result with a 
value of f pdv = 3/8 as in Bicknell et al. (1997), which we 
will adopt as a conservative lower limit. 

We numerically calculate the volume of the cocoon as 
it grows into the ICM and the pressure inside the cocoon at 
each step, and add up the pdV work to get the final amount 
of energy lost in this mode. The evolution of the cocoon is 
continued until the pressure inside the cocoon becomes equal 
to the ambient pressure, n a fcflT a (n a = p a /fj,m p ), where T a 
is the temperature of the ICM. The fraction f p dv is calcu- 
lated by taking the ratio of the total energy lost through 
mechanical work to the total energy (that is Ljtuf e ). (We 
found that for the relevant values of the ambient medium 
parameters, the time scale to reach pressure equilibrium is 
always larger than tuf e .) The dependence of f pdv on the 
ambient density, is shown in Figure 2 for T a = 10 6 , 10 7 and 
10 8 K, and for L k = 10 46 (solid lines) and 10 47 erg/s (dashed 
lines) . 

The plot shows that the fraction f p dv is a function of 
temperature of the cluster and also the density of the ambi- 
ent medium. The general trend is that the fraction reduces 
at higher temperatures and higher densities. This is because 
of the fact that the cocoon reaches pressure equilibrium with 
the ambient medium sooner for a higher pressure environ- 
ment (higher T a and n a ), and the total pdV work ends up 
being smaller. The plot also shows that the fraction f p dv 
depends weakly on the jet luminosity for lower tempera- 
tures (T a < 10 7 K) , whereas there is a bit of a difference for 
T a ~ 10 s K. 

The fraction f pd v calculated above is somewhat larger 
than that has been used in the literature, for ambient 
medium with low pressure. The difference is mainly the re- 
sult of our inclusion of cocoon evolution even after the jet 
has switched off. Incidentally, Inoue & Sasaki (2001), while 
using a value of f p dv = 1/4, discussed the possibility that 
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this fraction could be larger in reality, because of its contin- 
ued evolution after the switching off of the jet (their §3.2). 

In our calculation for the total pdV work done by BAL 
and RLQ outflows, we will use the values of f p dv obtained 
above. As mentioned earlier, the energetics and geometry of 
BAL outflows are not clear at present. For concreteness, we 
have worked out the case of RLQ outflows in detail, and we 
will use the same values of f p dv for BAL outflows as well. 



4 HEATING OF THE ICM 

Equipped with the knowledge of the rate of formation of 
quasars in clusters (equation O) and the fraction of total en- 
ergy which is deposited as pdV work by the outflows from 
them (§3.2), we are now in a position to calculate the total 
amount of non-gravitational energy provided by quasar out- 
flows in a cluster. If we denote the gas fraction of the total 
cluster mass by f ga s, then the total number of gas particles 
is ~ M c ifg aa /mp, where m p is the proton mass. It is not 
yet clear if the gas fraction has a universal value for clus- 
ters of all masses and at all redshifts. A few authors (see, 
e.g., Schindler (1999)) have found no correlation of the gas 
fraction with the cluster mass, whereas others (e.g., Ettori 
& Fabian (1999)) find some correlation (with low mass clus- 
ters possessing lower gas fraction) . This correlation has been 
attributed to the excess energy deposition (since low mass 
clusters are more liable to lose gas from heating than mas- 
sive clusters) (see, e.g., Bialek et al. (2000)). Here, we are, 
however, trying to calculate the magnitude of this very ex- 
cess energy. It would not be appropriate to include this cor- 
relation it a priori in our calculation. (We show later that 
including these correlations only increases our estimate of 
excess energy input.) We have, therefore, used a value of 
fgas =0.1 for our calculation. The total energy per (gas) 
particle deposited into the ICM of a cluster of mass M c i is 
then given by, 



tpdV 



where 



Mdfgas 

M cl 



u df q , c ,(M,z\M c i,zf) 
dMdz 



-j^-dM dz f [L k tufefpdv(n a ,T a )] . 



(17) 



d n q (M,z\M cl ,Zf) 



dMdz 



is calculated using equation 



|9|. The factor f a ~ 1 for the BAL outflows, and f ~ 0.1 for 
outflows from RLQs (§3). The redshift z m is the maximum 
redshift of heat input. We later show (Figure 7) that the heat 
input is negligible for z > 5. The density and temperature 
of the ICM of a cluster of a given mass (M c i) and formation 
redshift (zf) is calculated using (Eke et al. 1998), 



T a = 1.65xl0 7 (l+z) 
and, 

M c lfg as 



Mci 



W^h^Mp 



2 / 3 rfi A(fi ,; 



n{z) 



1/3 



where 



A(«) = 18tt + 82a: - 39a; 



,(18) 



(19) 



(20) 



and x = Q(z) — 1 (Bryan & Norman 1998), where we use 
equation [12 to compute Q,(z). For r„,: r , we use, 




log 10 (T) (keV) 

Figure 4. Excess energy (in keV) from BAL outflows is shown as 
a function of the cluster virial temperature (keV) for clusters with 
Zf = 0. The solid line shows the result of our calculation using 
the density and temperature dependent f p dv an d the dotted line 
shows the results when f p dv = 3/8. 



AirA(z)p cr it(z) 



1/3 



(21) 



where we have used M v i r = M c i, and p cr u is the critical 
density of the universe. The densities used in the following 
calculations range between 10" 4 -10" 6 cm . 

The integral in equation ( J]~t| ) is evaluated using Mi = 
10 12 M Q and M u = 10 13 M Q , for z f = 0,0.5,1 for dif- 
ferent values of M c i. We present the results for the total 
non-gravitational energy input per particle as a function of 
cluster mass (or, equivalently, gas temperature) in Figure 4 
(for zj = 0). The solid curve shows the heat input calcu- 
lated using fpdv from §3.2. The dotted line shows the case 
for a constant f p d v = 3/8 (Bicknell et al. 1997). We show the 
results for different Zf in Figure 5 (against T) and Figure 6 
(against cluster mass). 

We also show in Figure 7 the rate of deposition of en- 
ergy ( — j' ir ) as functions of the redshift for three clusters 
of masses M = 2 x 10 13 , 10 14 and 10 15 M , all for z f = 0. 
All the curves drop to zero at low redshift because of the 
condition dn q>c t/dt = ( in eqn ^). In reality the contribu- 
tion to the heating should be small but non-zero, and will 
increase the estimate of excess energy. 



5 DISCUSSION 

It has been estimated that the amount of excess energy re- 
quired to explain the observation is of order 0.5-3 keV per 
particle (Navarro et al. 1995; Cavaliere et al. 1997; Balogh et 
al.1999; Wu et al. 2000). Recently, however, Lloyd-Davies et 
al. (2000) have shown from observations of groups of clusters 
that an excess energy of 0.44±0.3 keV per particle suffices to 
explain the excess entropy in groups. They showed that this 
can explain the entropy floor for galaxy groups with tem- 
perature T ^ 4 keV. Borgani et al. (2001) have also shown 
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Figure 5. The excess energy is shown for Zf = (solid) 0.5 
(dotted) and 1 (dashed). 
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Figure 6. The excess energy against the cluster/group mass is 
shown for Zf = (solid) 0.5 (dotted) and 1 (dashed). 

with the aid of numerical simulations that excess energy of 
order ~ 1 keV per particle reproduces the observations. 

The solid and dashed curves in Figures 4 & 5 show that 
the excess energy from pdV work done by quasar outflows 
fall in this required range. It is seen that the excess en- 
ergy per particle is larger for clusters or groups with lower 
temperature. This is due to two factors: (a) the number of 
quasars per unit mass is larger for smaller clusters, and (b) 
the fraction of total energy in outflows that is lost in pdV 
work is larger for them. We show the results in the case of 
a constant f p dv = 3/8 (as in Bicknell et al. (1997)) with 
dotted lines. It is interesting to note that even in this case 
the excess energy is in the required range (0.44 ± 0.3 keV 



per particle), especially for groups with low temperatures, as 
advocated by Lloyd-Davies et al. (2000). Incidentally, this is 
larger than the estimate of excess energy from galactic winds 
(<, 0.1 keV per particle, Wu et al. (2000)). 

We found that our results for the excess energy (solid 
line in Figure 4) can be approximated by a fit of type (in 
keV per particle), 

/ T \ -0.193 / rr, \ -1.2 

^ ~°Kl0kev) +a033 (2kev) ' 

0.50 < T < 8.0 keV (22) 

We also found that the results for the excess energy taking 
fpdv =3/8 (Bicknell et al. (1997)) (dotted line in Figure 4) 
can be approximated by a fit of type (in keV per particle), 

/ T \ _0 ' 95 

e pdv ~ 0.17 + 0.0451 -— ,0.50 < T < 8.0 keV (23) 

V 2 ke V / 

We compare our results with the data from Lloyd- 
Davies et al. (200) in Figure 8. We show the predictions of 
our calculation as the thin solid curve (corresponding to the 
solid curve in Figure 4), where the data points have been 
taken from that of Lloyd-Davies et al. (2000) (their Figure 
9). We also show the result of the calculations for f p dv =3/8 
by the thick solid line. The data points refer to the binding 
energy of a (constant) central fraction (0.004) of the virial 
mass of groups and clusters. The dashed line shows the case 
for self-similar models, of type E oc T, derived from the data 
points for rich clusters. The dotted line shows their fit (with 
a constant excess energy of 0.44 keV per particle) along with 
a formal 1 a confidence interval shown by the shaded region. 
The figure shows that our predictions are consistent with the 
data available at present. The thick line (corresponding to 
fpdv = 3/8) falls close to the fit provided by Lloyd-Davies 
ct al. (2000), whereas the thin line (using f p dv from Figure 
3) somewhat overestimates the heat input at the low mass 
end. The thick line can be viewed as a conservative estimate 
of the heat input, since it uses f p dv = 3/8. The thin line, 
however, provides an estimate of the heat input if f pd v is 
much larger than 3/8. We should remind ourselves here that 
we have calculated f p dv for radio galaxies and used the same 
values for BAL outflows. If a more accurate estimate of f p dv 
for BAL outflows is worked out in the future, the resulting 
heat input into the ICM could then be scaled accordingly 
using Figure 8. 

We would like to emphasize here again that the our 
calculation provides a conservative estimate of the excess 
energy, for reasons outlined in §2. Moreover, we have used a 
constant density and temperature in time for the ICM gas 
(for clusters with a given z/), which is not very realistic. In 
reality, the density at higher redshift will be smaller, and 
the inclusion of a density dependent f p dv will only increase 
the estimate of excess energy (since this fraction increases 
with decreasing density). 

The curves for clusters forming at different epochs show 
excess energies to decrease somewhat for clusters with higher 
formation redshift. 

The curves of Figures 4 & 5 assume f a ~ 1, which is 
relevant for BAL outflows. The excess energy from RLQ out- 
flows will be one tenth of these curves, showing the difficulty 
of using radio galaxies as the only source of non-gravitational 
heating, if conservative estimates for their kinetic luminosi- 
ties are used. Recently, Inoue & Sasaki (2001) have used the 
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Figure 7. The rate of deposition of excess energy (—de p dv/d.z) 
is shown as a function of redshift, for M c ; = 2 X 10 13 Mq (solid 
line), 10 14 Mq (dotted line) and 10 15 M (dashed line) 

radio luminosity functions of Willot et al. (2001) and Ledlow 
& Owen (1996) to determine the abundance of radio galax- 
ies in clusters, and finally to estimate the total pdV work 
done by the cocoons of these radio galaxies. They estimated 
an excess energy of order 1 keV per particle for rich clusters 
like the Coma cluster, and also for poor groups, assuming 
that their ratio of radio galaxies per unit cluster mass is 
universal. From our calculation, we find an excess energy 
from only radio-loud quasars that is an order of magnitude 
smaller than their estimate. It is possible that the assump- 
tions leading to the estimate of Lu are at the source of this 
difference (see, e.g., the discussion on the uncertainty in the 
factor fj in their §3.2). 

The evolution of the rate of energy input (Figure 7) 
shows that the heat input before z ~ 5 is almost negligi- 
ble. Most of the heating occurs in the range z ~ 1-4. It 
is also clearly seen that the ICM of poor groups is heated 
at lower redshifts, compared to the gas in massive clusters. 
This follows from the simple consideration that the evo- 
lution of objects in a given mass range (here, that of the 
quasars) typically occurs earlier in more massive clusters. 
Our results, therefore, suggests that the ICM in clusters 
were "preheated", before the major mergers took place in 
them, whereas, the ICM in groups of galaxies were heated 
at epochs similar to that of their formation. It is interesting 
to note that the energy input epoch lies below the upper 
limit from recent observations of Lloyd-Davies et al. (2000). 
It is also consistent with the limit on the range of redshifts 
as shown by Yamada & Fujita (2001) given the level of un- 
certainty. 

Recently, accelerated particles from shocks as a result of 
the formation of clusters have been hypothesized to be the 
source for the diffuse gamma ray background (Loeb & Wax- 
man 2000). This gamma ray production will be suppressed, 
however, if the gas in the clusters were heated substantially 
at earlier epochs (Totani & Inoue 2001). Our results show 
that the ICM in massive clusters were preheated (and earlier 



than the ICM in groups), and the suppression of the gamma 
ray production will, therefore, be an important effect, if con- 
firmed. 

Finally, we discuss the uncertainties involved in our cal- 
culation. Apart from the uncertainties in cosmological pa- 
rameters, the major uncertainties lies in the factors f q (§2) 
(connecting the abundances of quasars with PS mass func- 
tion) , f (the fraction of quasars with outflows) ,f p dv and 
fgas- Among these, the most uncertain factor is f , which we 
have assumed to be of the order of unity for BAL outflows. 
The uncertainty in this factor will be reflected in the uncer- 
tainty of the final heat input (with a direct proportionality; 
see eqn |l^) . The uncertainty in f p dv has been already dis- 
cussed earlier, and we found that even if f p dv is as low as 3/8 
for all cases, the final excess energy is certainly larger than 
that from supernovae driven winds, and is still within the 
required range of excess energy, especially for loose groups. 
As far as the uncertainty in f gas is concerned, we have also 
done our calculation with a varying f ga3 , e.g., of the type, 

f gas = 0.15(1 + z)-°- 5 (Md/lO 15 !!- 1 ) 01 , (24) 

as has been advocated by Ettori & Fabian (1999), and we 
have found that the excess energy is approximately doubled 
in this case. It is, however, not clear if this correlation is a 
result of the excess energy, and, so, it would not be appro- 
priate to attach much significance to this result. We have 
also varied the lower limit in our estimate of f q (eq. [l]) and 
found that changing the lower limit from 10 12 Mq to 10 11 
M© increases the final heat input by only ~ 10 %. This 
is because of the fact that the increase in the number of 
quasars is compensated by the decrease in their mechani- 
cal luminosity. Lastly, we have already discussed in detail 
the uncertainty in the net formation rate of quasars in clus- 
ters, and as explained in §2, our approach here has been 
very conservative, and the final results should be regarded 
as conservative estimates in this regard. 



6 SUMMARY 

We have calculated the excess energy deposited by quasar 
outflows in clusters in order to explain the observations of 
excess entropy in groups and clusters of galaxies. We sum- 
marise our findings below: 

(1) We have used the extended Press-Sechter formal- 
ism to derive a formation rate of quasars inside clusters and 
groups, as a function of the cluster/group mass and its for- 
mation redshift. 

(2) We have calculated the fraction of the kinetic lu- 
minosity of outflows (RLQ and BAL outflows) that is de- 
posited onto the ambient medium, as a function of the den- 
sity and temperature of the ambient medium. For outflows 
from radio-loud quasars, we have included the evolution of 
the cocoon after the jet turns off. 

(3) The final excess energy from the mechanical work 
done by quasar outflows is found to be of order 0.18-0.85 
keV per particle, and is consistent with the data available 
at present. The excess energy in this scenario comes mainly 
from BAL outflows, with radio galaxies supplying about a 
tenth of the total. Keeping in mind the uncertainties in the 
estimate of energetics and abundances of radio and BAL 
outflows, we conclude that both radio galaxies and BAL 
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Figure 8. The prediction from our calculations is presented in the form of the final binding energy per 
particle of the central region of groups and clusters against the gas temperature. The data points are from 
the Figure 9 of Lloyd-Davics ct al. (2000), and the dashed line refers to their fit E oc T, derived from the 
data points for clusters with T > 4 keV. The dotted line refers to their second fit, with a constant excess 
energy of 0.44 keV per particle (subtracted from the binding energy)along with a formal 1 a confidence 
interval shown by the shaded region. The thick solid line uses f p dv = 3/8 and the thin solid line uses f p dv 
from Figure 3. 



outflows are promising candidates for heating the ICM. We 
found that this excess energy increases with decreasing mass 
of the cluster/group. This prediction could be tested with 
better data in the near future. The excess energy does not 
depend strongly on the formation redshift. 

(4) The epoch of heating is found to be in the range 
z ~ 1-4, where this epoch is at lower redshifts for low mass 
clusters. 
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